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Abstract 

Given a uniform ensemble of quantum density matrices p, it is useful to calculate 
the mean value over this ensemble of a product of entries of p. We show how to 
calculate such moments in this paper. The answer involves well known results from 
Group Representation Theory and Random Matrix Theory. This quantum problem 
has a well known classical counterpart: given a uniform ensemble of probability dis- 
tributions P = (Pi, P2, . . . , Pn) where the Pj are non-negative reals that sum to one, 
calculate the mean value over this probability simplex of products of P components. 
The answer to the classical problem follows from an integral formula due to Dirichlet. 
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1 Introduction 



The probability simplex of dimension N—l is the set of all points P = (Pi, P 2 , . . . , Pn), 
where the components of P are non-negative reals that sum to one. Take N = 3 for 
simplicity. It is useful to calculate moment integrals such as: 

J d 3 P P?P 3 , (1) 

where d 3 P = dPidP 2 dP3, and where the point P ranges over the 2 dimensional prob- 
ability simplex. As will be discussed in detail later, such integrals can be performed 
using an integral formula due to Dirichlet fl|] . 

Now suppose we generalize this problem to the quantum regime by consider- 
ing quantum density matrices instead of classical probability distributions. Suppose 
that p is a quantum density matrix (i.e., p is a Hermitian matrix with non- negative 
eigenvalues that sum to one). It is useful to calculate moment integrals such as 

J vp (ai) 2 pi, 2 • (2) 

As will be discussed in detail later, there is a very natural way of defining the measure 
Vp. The goal of this paper is to show how to calculate integrals like Eq.@. The 
answer involves well known results from Group Representation Theory @ {| f| ||, 
and Random Matrix Theory [[F] [|J [[|. Ref.|| by Itzykson and Zuber is especially 
pertinent to this paper. 

2 Notation 

In this section, we will introduce some notation that will be used in subsequent 
sections. 

RHS (ditto, LHS) will mean "right hand side" (ditto, "left hand side"). For 
any complex number z, let (ditto, Zq) represent its real (ditto, imaginary) part. 

As usual, 5(x) for real x will denote the Dirac delta function S(x) = Jf^ |^e ifcr . 
Likewise, 5(x, y) and b v x will denote the Kronecker delta function. b v x equals 1 if x = y 
and it equals if x ^ y. e^ia,...,^ will denote the totally anti-symmetric tensor with N 
indices, e^^...^ equals 1 if (z'i, i 2 , ■ ■ ■ , ijv) is an even permutation of (1, 2, 3, . . . , N), 
and it equals —1 for odd permutations. Let 9(S) be the "truth function" or "indicator 
function"; it equals 1 if the statement S is true, and it equals if S is false. For 
example, 9(x > 0) is the unit-step function; it equals 1 if x > and it equals if 
x < 0. 

The set of Hermitian NxN matrices will be denoted by Herm(N) . Any square 
matrix A is said to be positive semi-definite (ditto, positive definite) if the eigenvalues 
of A are non- negative (ditto, strictly positive). We will write A > (ditto, A > ) if 
A is positive semi-definite (ditto, positive definite). Suppose A and B are two NxN 



2 



matrices. We define the dot product of A and B by A ■ B = ti(A T B) = Mj^%,y 
We will use \ -J*— ]jj = ^ to denote the matrix of the partial derivatives with respect 



to the entries of A. Note that B ■ 



tr(B T ^ 



a A' 



We will often denote the set {xk ■ k 6 K} by {x k } keK , or simply by {x k }k 
when the set K is clear from the context. Likewise, the ordered set (x k '■ k e K) will 
be denoted by (xk)keK or simply by (x k )k- J2k&K x k will often be denoted by J2{ x k}k 
and UkeRXk by U{xk}k- 

Suppose a is an iV component column vector, a = (a , ai, . . . , «tv-i) T - One 
defines the Vandermonde determinant of a by 



A(o) 



«0 



On 



iV-l 



a 



N-l 



a 



AT-l 
iV-1 



(3) 



Often, we will express the RHS of the last equation more succinctly as follows: 

A(a) =det(a°,a\...,a 7V - 1 ) . (4) 
As is well known and easily proven: 

A(a) = II K ~ a i}o<i<j<N-i ■ ( 5 ) 

Vandermonde determinants are sometimes, especially in very old literature, referred 
to as difference products, or as basic alternants. 

We end this section by defining symbols for two expressions that arise fre- 
quently in subsequent sections. Let 



?n = 110 ! )i<i<7v • 
From Eq.(|), it follows that T N = A(0, 1,2, ... , N). Let 

N(N - 1) 



(6) 



(7) 



An N x N matrix has entries below (and above) its main diagonal. Also, 
l + 2 + -.. + (iV-2) + (JV-l). 



3 Ensemble of Classical Probability Distributions 

In this section, we will review the simplex moments integral and its variants. 

Suppose index b ranges over the integers from 1 to iV&. Define the following 
moment integral over an Nb — 1 dimensional simplex: 
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S Nb -t = l[{ r<ix b } ^x 6 -A)nW} 6 , (8a) 
where A is a positive real and the v\, are non-negative integers. It is well known that 

. um^ , (8b ) 

where 

v = £ v 6 + AT, - 1 . (8c) 

Eqs.(||) can be generalized so that they also apply to non- integer v h (just replace 
factorials by Gamma functions according to the prescription n\ — ► T(n + 1)), but we 
won't bother with such generalizations in this paper. Eqs.(^j) are a generalization to 
higher dimensions of the Beta function, given by 

ftm, n) = I' dt r-*(l - ty 1 = ^fr) . (9) 
jo 1 (m + nj 

Now let .B be an index that ranges over the integers from 1 to Nb, where 
Nb_ = iVft — 1. The iVg dimensional Dirichlet integral is defined by 

^-i=n{f° da *} c^< A )nw} fl /(E^)- ( ioa ) 

Uo Jb b r 



It is well known that 



where 



2Vi = , (10b) 



(7(A) = [ X dtf(ty- l v. (10c) 

Note that g(A) = A 1 ' when /(t) = 1. In Eqs.fllOD, v is defined as in Eq. (|8cD but with 
z/iVt = so that J2b v b = T,b v b- 

Proofs of Eqs. (||) for Sn b and Eqs. fllOD for £V S are easily found in the literature 
(see, for example Refs. 00) so we won't present them here. However, we do want 
to emphasize that these two integral formulas follow trivially from each other. One 
can prove as follows that the formula for V Nb implies the formula for Sn b - Take 4- 
of Eqs.([TDap and (gDTj ). Then use ±6(x > 0) = 5{x) and the fact that N R = N k - 1. 



Conversely, one can prove as follows that the formula for Sn b implies the formula for 

_ tnF,n fifth rri™ 



T> Nb . Applying the operator f2 = J A dX 4- to Eq.(^Fj) gives 



A ' dA /(A)%l = ni^f/ A 'dA/(A)r-0 , (li; 



dX v\ \Jo 
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whereas applying Q to Eq.(^a|) gives 

/;' dA/(A) £^i = (12a) 
= n{f * *}/ « /(A) ((-d4) *(i> - ^ <i2b) 

= n{^°°^BlB"} £ d\ f(\)6(£x B - \) (12c) 

= Dif t -i (12d) 

In going from line 6 to line c, we have assumed that z//v 6 = and performed the xjj b 
integration. 

4 Some Results from 

Group Representation Theory 

In this section, we will review quickly some well known facts from Group Repre- 
sentation Theory. For more details and proofs, see, for example, Refs. [§] @ 

B- 

A fT-box Young graph with at most N rows is specified by an iV-dimensional 
column vector of integers r\ = (rjo, r)i, . . . , rjN-i) T such that r/ > r/i > • • ■ > Vn-\ > 
and J2f=o l Vj = K- Vj is the number of boxes in the j-th row, where the top row is 
the 0-th one. Each subsequent row has the same number or fewer boxes than the row 
above. If some of the last few components of r\ are zero, they are often omitted. For 
example, 



in 



is specified by r\ = (6, 3, 3, 0, 0) T with K = 12 and N = 5. 

If Sjc is the symmetric group (permutation group) on K letters, then the 
classes of Sk are specified by a K-tup\e of non-negative integers i = (ii,i2, ■ ■ ■ 
such that lii + 2i 2 + . . . + Ki k = K. C(i) = (l n , 2* 2 , . . . , K %K ) represents the class of 
elements of Sk with i\ cycles of length 1, i 2 cycles of length 2, . . . , and %k cycles of 
length K. The order (i.e., number of elements) of C(i) is given by 

K\ 

' = (Pii 1 !)(2 j 22 2 !)...(^^!) ' ^ 13 ^ 

There is a one-to-one onto correspondence between: (l)the irreps of Sk, (2)the classes 
of Sk, (3)the Young graphs with K boxes. Tables |I] to £| give the characters for each 
(irrep, class) pair of Sk, where K ranges from 1 to 4. 

The irreps of the unitary group U(N) (and also the irreps of GL(N)) are in 
one-to-one onto correspondence with the Young graphs with at most iV rows. We 
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K=l 


(1) 


^class 


Characters\ 


1 


border 


□ 


1 




| irrep 







Table 1: Characters for each (irrep, class) pair of the permutation group Si. 



K 




=2 


(I 2 ) 


(2) 


^class 


Characters\, 


1 


1 


border 


m 


1 


1 










1 


-1 





| irrep 



Table 2: Characters for each (irrep, class) pair of the permutation group S 2 - 



will specify the irreps of U(N) by 77^, where 77 is an iV-dimensional column vector of 
integers that specifies a Young graph. (The superscript N serves to distinguish this 
from an irrep of a symmetric group). In terms of tensors, the number K of boxes of the 
Young graph corresponds to the number of tensor indices, {1, 2, ... , N} corresponds 
to the range of the tensor indices, and the Young graph gives the symmetry properties 
of the tensor. 

Weyl showed that the dimension of the irrep r] N of U (N) is given by: 

dim(^) = d<*«r, + 6) N -\(r, + *)*-*,..., (V + S)°) {u) 

where 5 = (N - 1, TV - 2, . . . , 1, 0) T . 

Let x v (A) represent the character of A e U(N) in the irrep r] N of U(N). 
Inspired by an identity due to Frobenius, Weyl derived the following two expressions 
for x vN (A). First, 

X 1 } det(a N -\a N - 2 ,...,a°) ' 1 j 

where a is the vector of eigenvalues of the N x N matrix A. Second, 

X"V) = E ^X v {C{m% 2 ■ ■ , (16) 
C(i) A • 

where K is the number of boxes in the Young graph 77, and where the sum is over all 
classes C(i) of the permutation group Sk, \ C(i) \ is the order of class C(i), x v {C(i)) is 
the character of C(i) in the irrep 77 of Sk, and where ti, t 2 , ■ ■ ■ are defined by 
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K= 


=3 


(I 3 ) 


(1,2) 


(3) 


^class 


Characters\ 


1 


3 


2 


border 


1 1 1 1 


1 


1 


1 










2 





-1 










1 


-1 


1 





t irrep 



Table 3: Characters for each (irrep, class) pair of the permutation group S3. 



K= 


=4 


(I 4 ) 


(I 2 , 2) 


(1,3) 


(2 2 ) 


(4) 


<— class 


Characters\, 


1 


6 


8 


3 


6 


<— order 





1 


1 


1 


1 


1 








1 1 


3 


1 





-1 


-1 










2 





-1 


2 













3 


-1 





-1 


1 










1 


-1 


1 


1 


-1 





I irrep 



Table 4: Characters for each (irrep, class) pair of the permutation group S4. 



N-l 

tr(ZT) = £ {a 3 Y , (17) 

j=0 



where {ctj}j are the eigenvalues of A. 

Eq.(14) for dim(r] N ) follows from Eq. flTo]) and the fact that dim(r/ JV ) = x v (I), 
where / is the identity matrix. Let RHS1 denote the RHS of Eq.(|15|) and RHS2 denote 
the RHS of Eq. (|i~6l) . Frobenius was the first to prove that RHS1=RHS2, but his proof, 
which is discussed in Hamermesh[[|], made no mention of U(N). Weyl gave a new 
proof in which U (N) was crucial. 

Note that in Eqs. QT5D and ([TB]) , A is a unitary matrix. Hence, its eigenvalues 
are of the form <x, = e lZ - Qj for some real Z<x,. However, the equation RHS1=RHS2 can 
be analytically continued to complex atj with \ctj\ 7^ 1. This is because both RHS1 
and RHS2 are linear combinations of monomials of the form a^af 1 . . . oi^Si , where 
PQiPii ■ ■ ■ jPn-i are non-negative integers. Since the Za/s are arbitrary reals, the 
coefficient in RHS1 of any fixed monomial a^a^ 1 . . . a^Si must equal the coefficient 
in RHS2 of that same monomial. 
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Table [5] was calculated using Eq.([T(J). A very similar table can be found in 
Ref. ||. Table ^was derived from the information in Table [|. 



K 


irrep 


X nN (A) 


dim( V N ) =x v "(I) 


1 


□ 


*1 


N 


2 


m 


lj-2 lj- 
2*i + 2*2 


±N(N + 1) 






1,2 1 4- 
2*1 ~~ 2 t2 


\N(N - 1) 


3 


i i i i 


1 J-3 l 1 -/- -/- 1 1 4- 

qZ 1 + 2*1*2 + 3*3 


%N(N + 1)(N + 2) 






\t~l + 0*i*2 — i*3 


iJV(iV + l)(JV-l) 






|*f — 1*1^2 + |*3 


|Ar(Ar-l)(Ar-2) 


4 





1412 12 1 1 
24*1 + 4*1*2 + 8*2 + 3*1*3 + 4*4 


^Ar(AT+l)(AT + 2)(Ar + 3) 




i i 


8*1 + 4*1*2 — 1*1 + 0*1*3 — 4*4 


lN{N + l)(N + 2){N-l) 






^*1 + 0*^*2 + ^*2 — |*1*3 + 0*4 


±N 2 (N + 1){N -1) 






8*1 — 4*1*2 ~ |*2 + 0*1*3 + i*4 


lN{N+l)(N-l)(N-2) 






2l*l _ 1*1*2 + |*2 + 1*1*^3 — |*4 


±N(N-l)(N-2)(N-3) 



Table 5: K is the number of boxes in the Young graph. The irreps of U(N) and GL(N) 
are in 1-1 correspondence with the Young graphs with at most A" rows. x v {A) is the 
character of A e U(N) in the irrep r] N of U(N). dim(r] N ) is the dimension of irrep 



5 Ensemble of Quantum Density Matrices 

In this section, we will generalize the simplex moments integral Eqs.(|8|) to the quan- 
tum realm. To go from classical to quantum physics, we will replace probability 
distributions by quantum density matrices. In the quantum case, we will need to 
do integrals over a manifold of matrices. Such integrals are used in several fields of 
mathematical physics. They are crucial to the field of Random Matrix Theory . 
Let 

It is easy to check that U is unitary. For any complex number z, let z^ (ditto, zq) 
represent its real (ditto, imaginary) part. Note that 
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K 


Sail irreps r) N with K boxes dim^^)^ (A) 





1 (convenient definition) 


1 




2 




3 


%(NH 3 1 + 3Nt l t 2 + 2t 3 ) 


4 


%(N 3 4 + 6N%t 2 + 8Nt 3 ti + 3Nt 2 + 6t 4 ) 



Table 6: This table was derived using the information in Table 



^ Ul ^ I • (19) 



Z I \ Zc± 



V2 



This transformation rule motivates us to generalize the volume element dx in real 
space to a volume element d 2 z in complex space according to: 

dx — > <i 2 z = 2dzsftdz<s . (20) 

We also generalize the Dirac delta function S(x) in real space to a Dirac delta function 
<5 2 (z) in complex space: 

S(x) - 5 2 (^) = ^(^)<5(z s ) . (21) 

Thus, 

= 1 -> d 2 zS 2 (z) = l. (22) 



Next we want to give a convenient definition of a volume element £>p for the 
manifold of all Hermitian N x N matrices p. We define 

V P = U{dp**hI[{d 2 P**'} X!X , x<xl ■ (23) 
This definition of T>p prompts us to define: 

s(p)=n{5(p~)},n{pM}^ tx<IBl , (24) 



so that 



/ Vp 8{p) = 1 . (25) 



The manifold of real points x with volume element dx corresponds to the 
manifold of Hermitian N x N matrices p with volume element T>p. But what is the 
matrix counterpart of the manifold of complex points z with volume element d 2 z7 
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A natural candidate for this is the manifold of all complex N x N matrices A. We 
define its volume element by: 



V 2 A = il{d*A xx ,} xxi . (26) 

Note that V 2 A is a product of twice as many dxASke real-space volume elements as 
T) 2 p. Finally, we define 

5\A)=\[{5\A xx ,)} xxi (27) 

so that 

J V 2 A 5 2 (A) = 1 . (28) 
Given two Hermitian N x N matrices p and u, one has 

UJ ■ p = ^LOxxPxx + 2 ^ [( U xv)si(Pxy)?ft - (^xy)^(Pxy)^} • (29) 
x x<y 

Combining the last equation and the identity 5(x) = Jf^ ^e lkx for real x, yields 

Vuj exp(zcj ■ p) = (2n) N2 5(p) . (30) 



Let X be a Hermitian N x N matrix with eigenvalues {Xj}j- We need to 
consider a real valued function f(X) that is invariant under unitary transformations 
of its argument X; that is, f(UXU') = f(X) for any unitary N x N matrix U. 
f(X) depends only on the eigenvalues of X, so f(X) = F(x) for some function 
F : H N — > R. Note that not all functions from to R qualify for the job of 
F. F(x) must also depend on x in a symmetrical way: since all N x N permutation 
matrices belong to U(N), F(x) must be invariant under permutations of its arguments 
{Xj}j- Henceforth, we will indulge in a convenient abuse of notation by replacing the 
symbol F by / so that f(X) = /(%). 

Let X and A be Hermitian N x N matrices with eigenvalues {Xj}j an d 
respectively. Let f(X) be a real valued function of X that is invariant under unitary 
transformations of its argument X. Finally, let £ be an arbitrary complex number. 
The following two integral formulas are well known: First (see Ref.[ 

(27TJ 

~N ~~ {J—oc ) j 



f vx f(x) = n { r d A A2 (x)m , (31) 

J J~N <J — oo ) j 



and, second(see Refs. 



/*«««^M=(^n{£<^tg^-vcx) ■ (32) 
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In the spirit of non-rigorous, applied mathematics, I will not specify precise sufficiency 
conditions on A, /(■) and £ under which Eqs.(|3~l~D and (|32| ) are valid. I leave it to 
more competent pure mathematicians to figure this out. 

Proofs of Eqs.(PTD and are easily found in the literature so we won't 
present them here. Proving Eq . (|32|) involves solving a partial differential equation 
(an initial value problem for a diffusion equation). 

Of course, Eq.(^TJ) follows from Eq.(p2|) when £ — > 0. To check this quickly, 
assume (without loss of generality) that the eigenvalues atj are very small and widely 
separated: 

< a < ax < ■ ■ • < a N _i < 1 . (33) 
Then, keep the largest term in the expansion of A(a): 

A(a) « a ai . . . a^-i • (34) 

Furthermore, keep only the largest term in the Taylor expansion of exp(£x • a) and 
in the multinomial expansion of (x ■ o) Ln : 



exp(£x • a) « — « — -— (xoXi • • • Xn-i)\OqOix ■ • -ajv-i) • (35) 

L N \ L N \J- N _i 

Use these approximations on the integrand on the RHS of Eq. fl32|) . Also replace the x 
dependent part of the integrand, that is, A(x)xoXi ' ' ' Xn-i, by its totally symmetric 
part A 2 (x)/N\. This converts the RHS of Eq.® into the RHS of Eq.(^) in the 
limit £ -> 0. 

We are almost ready to present our generalization of the simplex moments 
integral. But first we need to prove two lemmas. 

Lemma 5.1 Suppose /?o, • • • ,Pn-i ar ^ non-negative integers. Then 

II + = II (A !},• A(/3) (36) 

proof: 

Note that the LHS of Eq.([36|) can be expressed as det(M), where the matrix 
M has entries Mjj = (z + Eq.(^) can be proven easily using mathematical 

induction and simple properties of determinants. QED 

Lemma 5.2 Suppose Po, Pi, ■ ■ ■ , Pn-i ar e non-negative integers. Then 

n {/; ^ - . . . & = ( I 1} , (37) 
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proof: 

Replace A(x) by 

A(x) = 6fcoA;i---fcjv-i a '0 0a 'l 1 • • ' X N—1 • (3^) 

in the LHS of Eq.(^). Then use the simplex moments integral Eqs.(|J) to get 

lhs = ^■■■^^ + ^ l h m 

where 

v = +k j )+N-l = y £p j + L N + N-l. (40) 

j 3 



Finally, apply Eq.fl3|). QED 

The volume of a simplex [v^ = V6 in Eq.(§) ) is the simplest case of the 
simplex moments integral. To warm up, we first present the quantum version of this 
simplest case: 



Claim 5.1 Suppose X e Herm(N). Then 

J VX 9(X > 0)6(trX - 1) = ■ (41) 



proof: 

Let LHS (ditto, RHS) stand for the left (ditto, right) hand side of Eq.@- 

Then 



LHS = n { r d A A2 (^)^(E X 3 - 1) (42a) 

{27T)LN -N\ I] { jH d X] ) . Hx)x°oX\ ■ ■ ■ xJfcfrE X 3 - 1) (42b) 

(42c) 



•Fn 

(27r) L -A(0,l,2,...,iV-l) 



(2L N + N- 1)! 

= RHS. (42d) 

In going from line a to line b, we replaced one of the A(x) by NIxqXi ■ ■ ■ Xn-v this was 
valid because the rest of the integrand was totally ant i- symmetric under permutations 
of the {Xj}j- To go from line b to line c, we applied Eq. (|37|) . To go from line c to line 
d, we used A(0, 1, 2, . . . , AT — 1) = F N _ X . QED 

Finally, we are ready to present the main result of this paper, a generalization 
of the simplex moments integral to quantum mechanics. Actually, we will give a 
moment generating function and calculate moments from that. 
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Claim 5.2 Suppose X,Ae Herm(N) , and £ is a complex number. If both sides of 
the following formula exist, then 



J VXe^9(X > 0)*(trX-l) = £ ^ Z^-lY. E dim(^) X ^(A) , 

K=0 v /' rj N :K boxes 

(43) 

where we define dim(rj N )x v {A) = 1 /or = 0. 
proof: 

We will prove this claim when £ = 1. The more general case can be obtained 
from this by scaling A (i.e., replacing A by £A, where £ is real, and then analytically 
continuing £ to complex values.) Let LHS (ditto, RHS) stand for the left (ditto, right) 
hand side of Eq.©. 

Applying Eq.(^) yields 

LHS = (2tt)^ J] { /_~ . ^ye^ - J] {e( Xj > 0)}, 6(£ Xj ~ 1) • (44) 
Recall the multinomial expansion: 

N K\ 

! + x 2 + - ■ ■ + x N ) K = E 5(E k r, K)———x^x k ^ ■ ■ ■ x k N N . (45) 



[x 



k 



nw} T 



Using the Taylor expansion of exp(-) and the multinomial expansion, we get 



00 (y ■ a) K> 

exp( X -a) = £ ^— L_ (46a) 

if'=0 - 



EE^E^^nj^g^} ■ (46b) 



X'=0 % r \ 

Now apply to Eq.(f44j) this expansion of exp(x • a) and Eq. (p?[) : 



lhs = f E'S>, n R} r E^^tfU-D, ■ <«) 



Jf'=0 fe 



The function being summed over A; is a product of a totally anti-symmetric 
function of k times n j ^} ■ Hence we may replace n{ a r r } by its totally anti- 
symmetric part: 
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II "> det (« fc °' « fcl > • • • . « fcAr_1 ) • (48) 

After doing this we will have a sum over k of a totally symmetric function of so we 
can replace: 

£ _> jvi £ 0(0 < k < fa < ■ ■ ■ < k N -x) . (49) 

k k 

(Terms with k{ = ki+i vanish.) We can also change the lower limit of the K' sum 
from K' = to K' = Ln, because terms with K' < Ln do not contribute. Let 
K = K' — Ljy- Let us change variables from K' to K. All these changes yield: 



7«C 



LHS = (27r)^EE 

A'=0 I 

6(0 < k < fa < ■ ■ ■ < fcjv_i)(JE K + Ln ' 



A{k) det(a k °,a k \... 



a 



(K + N 2 -l)\ det(o°, a 1 , . . . , a"" 1 ) 
Next consider the following change of variables: 

?7V-2 = fa ~ 1 

771 = k N _ 2 — (iV — 2) 
770 = fci\r-i — (iV — 1) 



(50) 



(51) 



Note that rjo — rji = fcjv-i — fc/v-2 — 1 > 0, 771 — ?72 = — — 1 > 0, etc. Also, 
J2j t]j = K' — Ln = K. Hence, 77 specifies a Young graph with K boxes and at most 
iV rows. 

Define a column vector 5 = (N - 1, N - 2, • • • , 1, 0) T . Let R be the N x N 
matrix that has ones on the non-principal diagonal and zeros everywhere else. For 

example, for N = 2, R = ( ^ q ^ ^ c ^ mens * ona ^ c °l umn vector 77, define 

its R transform r] R by rj R = Rr]. R just reverses the entries of rj. Eqs. ([5l|) can now be 
stated succinctly as 77 = k R — 5. Note that R 2 = 1. Thus, det(RAR) = det(A), for any 
square matrix A of the same dimension as R. The right and left i?'s reverse the order 
of the columns and of the rows of A, but this does not change the value of det(A). For 
example, suppose det(A) is the Vandermonde determinant A(v) of a vector v. Assume 
v is 3 dimensional for concreteness. Then det(v°, v 1 , v 2 ) = det((v R ) 2 , (v 11 ) 1 , (v R )°). 
Changing the summation variable from k to 77 in Eq.(|50|) yields: 
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MS = (JV-l)l +(2 '> gE (if + iv2 _ 1)! 7, (52) 
where 7 is defined by 

_ det(a riN -\a m -2 +1 , . . . , a™^- 1 ) 
7 ~ detfa ,^ 1 ,...,^- 1 ) ' ( } 



Reversing the order of the columns of both the denominator and numerator 
determinants of 7 yields: 



7 = "^1^" ' = **" (A) . (54) 



det (a'» +N - 1 , a^+ N - 2 , a™-^ 
det{a N -\a N - 2 ,...,a° 

Reversing columns and rows in A(r] R + S R ) leads to Weyl's formula for the dimension 
of 77^. 



A( V R + 5 R ) = det[( v + 5) N ~\ ( V + 5) N ~ 2 , . . . , (77 + 5)°] = dim( v N )F N ^ . (55) 



Applying Eqs.fl54|) and (|55]) to Eq.(||) finally yields LHS= RHS. QED 

Eq.(fi3|) gives a generating function that can be used to calculate moments 
over a uniform ensemble of density matrices. For example, we can calculate the 
mean value of X ix ^X i2 ^ 2 . Such moments will have indices attached because X is a 
matrix. To avoid having free indices, we will bind them to constant matrices. So 
instead of calculating the mean value of Xi 1 j 1 Xi 2 j 2 , we will calculate the mean value 
of (Ci • X)(C*2 • X), where C\ and C 2 are constant matrices. 

For any N x N matrices A and C, consider how the operator C ■ ^ acts on a 
power A n for some integer n. C • ^ replaces one A at a time by a C: 

d n ~ l 

C ■ —A n = y A'CA"- 1 -* . (56) 



For example, if we define the operator Q by 



then: 

fit* = ^tr(C 1 )tr(C 2 )tr(C 3 )tr(C 4 ) , (58a) 
p 

nt 2 tj = ^tr(C 1 C 2 )tr(C 3 )tr(C 4 ) , (58b) 
p 

m==X;tr(C 1 C7 2 )tr(C3C4), (58c) 
p 
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fit 3 ti = E tr(C 1 C 2 C 3 )tr(C , 4 ) , (58d) 
p 

Qt 4 = Etr(CiC 2 C3C 4 ) . (58e) 
p 

In Eqs.fl5"8"|), t r . = tr(v4 r ) as before, and the sums run over all permutations P on 4 
letters. P acts on the subscripts {1,2,3,4}. 
Define 

Z{A) = J VX e tl{AX) 6(X > 0)5(trX - 1) , (59) 

1(d) = JvX (d ■ X)6(X > 0)8(trX - 1) , (60) 

I(d, d) = JVX [d ■ X)(d ■ X)6{X > 0)5(trX - 1) . (61) 
Then, by virtue of Table |] and Eg. (^3|) , one has 

1(d) = lim (d ■ 4j) Z ( A ) ( 62a ) 



a^o\ OA 
{2<k) l »F n _ x 



(N 2 )\ 

and 



iVtr(d) , (62b) 



I(d,d) = ^.^(d-^ZiA) (63a) 



{27t) Ln ^ 



"^NlNtridMd) + tr(dC 2 )} ■ (63b) 



(N 2 + 1)! 
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